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ABSTRACT 

Understanding the origin of stellar masses is a key problem in astrophysics. In the solar 
neighborhood, the mass distribution of stars follows a seemingly universal pattern. In 
the centre of the Milky Way, however, there are indications for strong deviations and 
the same may be true for the nuclei of distant starburst galaxies. Here we present the 
first numerical hydrodynamical calculations of stars formed in a molecular region with 
chemical and thermodynamic properties similar to those of warm and dusty circum- 
nuclear starburst regions. The resulting IMF is top-heavy with a peak at ~ 15 M , 
a sharp turn-down below ~ 7 M Q and a power-law decline at high masses. We find 
a natural explanation for our results in terms of the temperature dependence of the 
Jeans mass, with collapse occuring at a temperature of ~ 100 K and an H2 density of a 
few xlO 5 cm" 3 , and discuss possible implications for galaxy formation and evolution. 
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1 INTRODUCTION 

Identifying the physical processes that determine the masses 
of stars and their statistical distribution, the initial mass 
function (IMF), is a fundamental problem in star-formation 
research. It is central to much of modern astrophysics, with 
implications ranging from cosmic re-ionisation and the for- 
mation of the first galaxies, over the evolution and structure 
of our own Milky Way, down to the build-up of planets and 
planetary systems. 

Near the Sun the number density of stars as a func- 
tion of mass has a peak at a characteristic stellar mass of a 
few tenths of a solar mass, below which it declines steeply, 
and for masses above one solar mass it follows a power-law 
with an exponent dN/d\ogm oc m -1,3 . Within a radius of 
several kpc this distribution shows surprisingly little varia- 
tion (Salpeter 1955; Scab 1998; Kroupa 2001; Kroupa 2002; 
Chabrier 2003). This has prompted the suggestion that the 
distribution of stellar masses at birth is a truly universal 
function, which often is referred to as the Salpeter IMF, al- 
though note that the original Salpeter (1955) estimate was 
a pure power-law fit without characteristic mass scale. 

On the other hand, there is increasing evidence that the 



IMF close to the centre of our Milky Way (Stolte et al. 2002, 
2005, Nayakshin & Sunyaev 2005, Paumard et al. 2006) and 
the neighboring Andromeda galaxy (Bender et al. 2005) is 
dominated by massive rather than low-mass stars. For the 
circum-nuclear starburst regions in more distant galaxies, 
very similar IMF deviations are subject to continuing debate 
(e.g., Scalo 1990, Elmegreen 2005). However, no conclusion 
has yet been reached, and it appears timely to examine the 
problem from a theoretical point of view. 

We approach the problem by means of self-consistent 
hydrodynamical calculations of fragmentation and star for- 
mation in interstellar gas where chemical and thermody- 
namical properties are described by a realistic equation of 
state (EOS). We focus on the most extreme environmental 
conditions such as occur in the nuclear regions of massive 
star-forming spiral galaxies. There the inferred dust and gas 
temperatures, gas densities and star formation rates typi- 
cally exceed the solar-neighborhood values by factors of 3, 
10 and > 100, respectively (e.g. Ott et al. 2005; Israel 2005; 
Aalto et al. 2002; Spinoglio et al. 2002). Consequently, it has 
long been speculated that such conditions lead to deviations 
from the Salpeter IMF (e.g., Scalo 1990, Elmegreen 2005). 
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Figure 1. Column density distribution of the gas projected along the three principal axes of the system after 3.5 million years of evolution 
when 15% of the gas is converted into protostars (sink particles). Their location is indicated by black dots. 



2 MODEL 

Stars and star clusters form through the interplay between 
self-gravity on the one side and turbulence, magnetic fields, 
and thermal pressure on the other (for recent reviews see 
Larson 2003; Mac Low & Klessen 2004; Ballesteros-Paredes 
et al. 2006). The supersonic turbulence ubiquitously ob- 
served in interstellar gas clouds can create strong density 
fluctuations with gravity taking over in the densest and most 
massive regions. Collapse sets in to build up stars and star 
clusters. Turbulence plays a dual role. On global scales it 
provides support, on local scales it provokes collapse. Stel- 
lar birth is thus intimately linked to the dynamic behavior 
of the parental gas cloud, which governs when and where 
star formation sets in (as illustrated in Figure 0. 

The chemical and thermodynamic properties of inter- 
stellar clouds play a key role in this process. In particular, 
the value of the polytropic exponent 7, when adopting an 
EOS of the form P oc p 7 , strongly influences the compress- 
ibility of density condensations as well as the temperature 
of the gas. The EOS thus determines the amount of clump 
fragmentation, and so directly relates to the IMF (Vazquez- 
Semadeni et al. 1996) with values of 7 larger than unity lead- 
ing to little fragmentation and high mass cores (Li, Klessen, 
& Mac Low 2003; Jappsen et al. 2005). The stiffness of the 
EOS in turn depends strongly on the ambient metallicity, 
density and infrared background radiation field produced by 
warm dust grains. The EOS thus varies considerably in dif- 
ferent galactic environments (see Spaans & Silk 2000, 2005 
for a detailed account). 

For the circum-nuclear starburst regions that are the 
subject here we assume a cosmic ray ionisation rate of 

3 x 10~ 15 s -1 , solar relative abundances (Asplund et al. 
2004; Jenkins 2004) and an overall metallicity of two times 
solar (Barthel 2005). A velocity dispersion AKur = 5 km/s 
is adopted to take the larger input of kinetic energy (e.g. 
through supernovae) into account. The dust temperature in- 
side the model clouds is set by a fiducial background star for- 
mation rate of 100 Mq yr _1 /kpc 2 which causes dust grains 
to be at temperatures of about Td — 30 — 90 K, depending 
on the amount of shielding. Gas temperatures range from 
T g = 40 - 140 K, over a density range of 10 4 - 10 7 cm~ 3 . 
These values are consistent with gas and dust temperatures 
determined for circum-nuclear starburst regions (Klaas et 
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Figure 2. Starburst EOS adopted from Spaans & Silk (2005). 



al. 1997; Aalto et al. 2002; Spinoglio et al. 2002; Ott et al. 
2005; Israel 2005). 

Figure [5] shows the resulting polytropic exponent as a 
function of density. The main feature is the 7 > 1 peak 
around n = 10 4 cm 3 . This peak implies that the gas warms 
up as it is compressed and it is caused mainly by strong pho- 
ton trapping in opaque H2O and CO lines in the metal-rich 
nuclear gas. That is, the large optical depth in the cool- 
ing lines suppresses the cooling efficiency. Also, warm dust 
(T > 40 K), heated by the ambient stars, causes H2O col- 
lisional de-excitation heating through far-infrared pumping 
(Takahashi, Hollenbach & Silk, 1983; Spaans & Silk 2005), 
which adds to the gas-dust heating. Cosmic-ray heating rate 
is elevated by a high supernova rate, as expected for nuclear 
starburst regions (Bradford et al. 2003). 

Adopting this EOS we follow the dynamical evolution 
of the star-forming gas using smoothed particle hydrody- 
namics (SPH). This is a Lagrangian method to solve the 
equations of hydrodynamics, where the fluid is represented 
by an ensemble of particles, and flow quantities are obtained 
by averaging over an appropriate subset of SPH particles 
(Monaghan 2005). The method is able to resolve high den- 
sity contrasts as particles are free to move, and so the parti- 
cle concentration increases naturally in high-density regions. 
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The performance and convergence properties of SPH are well 
understood and tested against analytic models and other nu- 
merical schemes in the context of astrophysical flows (see, 
e.g., Mac Low et al. 1998; Lombardi et al. 1999; Klessen et 
al. 2000; O'Shea et al. 2005; Ballesteros-Paredes et al. 2006). 
Artificial fragmentation can be ruled out, as long as the mass 
within one smoothing volume remains less than half the crit- 
ical mass for gravitational collapse (Bate & Burkert 1997; 
Hubber, Goodwin, & Whitworth 2005). We use the publi- 
cally available parallel code GADGET (Springel et al. 2001). 
It is modified to replace high-density cores with sink parti- 
cles (Bate, Bonnell, & Price 1995) that can accrete gas from 
their surroundings while keeping track of mass and momen- 
tum. This enables us to follow the dynamic evolution of the 
system over many local free-fall timescales. We identify sink 
particles as the direct progenitors of individual stars. For a 
more detailed account of the method and a discussion of its 
convergence properties we refer the reader to Klessen et al. 
(2000) and Jappsen et al. (2005). 

We focus on a cubic volume of 11.2 pc in size, which 
contains 80, 000 Mq of gas and has an initial mean particle 
density n = 10 3 cm 3 at a temperature of 21 K. Above the 
characteristic density n = 10 4 cm -3 where 7 is at a max- 
imum, the temperature quickly reaches values of ~ 100 K. 
This set-up is chosen to describe the typical environment 
within the central regions of an actively star-forming galaxy 
such as our own Milky Way or NGC 253. In such galaxies, 
high-density gas with n > 10 5 cm , as traced by HCN, 
typically has a filamentary structure with very low filling 
factor, while the bulk of the gas is at n ~ 10 3 cm -3 (Mor- 
ris & Serabyn 1996; Hiittemeister et al. 1993; Israel & Baas 
2003), exactly as found at the end of our simulation (see 
Figure 0. We stop the calculation at a star formation effi- 
ciency SFE « 15%, when roughly 1/6 of the total gas mass 
has turned into gravitationally collapsed condensations (i.e. 
sink particles, which we identify as direct progenitors of in- 
dividual young stars). Throughout the simulation we drive 
turbulence continuously on large scales, with wave numbers 
k in the range 1 < k < 2 (see Mac Low 1999) to yield a 
constant turbulent Mach number 7Vf rms ~ 5. The particle 
number is N — 8 000 000. This is thus one of the highest- 
resolution star-formation calculations done with SPH, with 
a total CPU time of 8 x 10 4 hours. The critical density for 
sink particle formation is n c = 10 7 cm~ 3 , with a sink parti- 
cle radius of 0.015 pc. The mass of individual SPH particles 
is m = 0.01 Mq, which is sufficient to resolve the minimum 
Jeans mass in the system Mj ~ 1.5 Mq. Except for the EOS 
and the particle number, the numerical set-up is identical to 
the study by Jappsen et al. (2005). We have performed a 
second run for a region of 5.7 pc with 4 times less mass, 
eight times fewer particles and a sink particle radius of 0.02 
pc that has reached a SFE ~ 36%. 



3 RESULT AND PHYSICAL 
INTERPRETATION 

We find that in the considered star-forming region, the mass 
spectrum of collapsed objects is biased towards high masses. 
The resulting IMF has a broad peak at ~ 15 Mq followed by 
an approximate power-law fall-off with a slope in the range 
-1.0 to -1.3. Furthermore, there is a clear deficit of stars be- 



low 7Mq. This is illustrated in Figure^. We contrast this 
finding with the result from a simulation appropriate for 
the physical conditions in star forming regions near the Sun 
(from Jappsen et al. 2005), where 7 changes from 0.7 to 1.1 
at an H2 density of a few xlO cm" . As expected, Figure 
shows a mass spectrum that is very similar to the IMF 
in the solar neighborhood (Kroupa 2002; Chabrier 2003). 
These striking differences are caused by the very disparate 
chemical and thermodynamic state of the star forming gas 
in the two simulations, since all other parameters are very 
similar. Our results thus support the hypothesis that for ex- 
treme environmental conditions as inferred for the centres of 
most spiral galaxies or more general for IR-luminous circum- 
nuclear starburst regions the IMF is indeed expected to be 
top-heavy. 

There is a natural explanation for our results in terms 
of the temperature dependence of the Jeans mass Mj . Com- 
pared to a mean temperature of 10 K for dense molecular 
gas in the Milky Way, gravitationally collapsing gas in our 
simulations has a temperature of ~ 100 K and an H2 density 
of a few xlO 5 cm -3 . As the critical mass for gravitational 
collapse scales as Mj oc T 1 ' 5 , this boosts Mj from 0.3 Mq 
at 10 K to about 10 M at 100 K (see also Klessen & Burk- 
ert 2000, Bonnell, Clarke, & Bate 2006). This temperature 
may seem high, but is quite consistent with molecular cloud 
observations in the Galactic centre (e.g. Hiittemeister et al. 
1993) or with high-density (n > 10 4 cm" 3 ) NH 3 data in the 
starburst centre of NGC253 (Ott et al. 2005). We also note, 
that this Jeans mass scaling argument is supported by recent 
observations in more nearby high-mass star-forming regions. 
For example, in M17 at a distance of 1.6 kpc from the Sun, 
the mass spectrum of prestellar cores, which are the direct 
progenitors of individual stars, peaks at at ~ 4Mq at an 
ambient temperature of 30 K (Reid & Wilson 2006) . This is 
well above the corresponding peak in low-mass star-forming 
regions (e.g. Motte, Andre, & Neri 1998). 



4 DISCUSSION 

Our mass spectrum is in good agreement with the IMF es- 
timates in the Galactic centre by Stolte et al. (2002, 2005), 
Nayakshin et al. (2005), and Paumard et al. (2006). For ex- 
ample, Stolte et al. (2002, 2005) find for the Arches cluster 
a clear deficit of stars below 7Mq. This is consistent with 
our result in the sense that the ambient densities and tem- 
peratures found in the Galactic centre are similarly elevated 
(Heifer & Blitz 1996) as in the circum-nuclear starburst en- 
vironment we consider. We stress that the turn-down in our 
model IMF at masses below 10 Mq is a direct consequence of 
the stiff EOS for densities n above a fewxlO 3 cm" 3 through 
the Jeans mass temperature dependence, and is not caused 
by resolution effects. Our two simulations resolve masses 
down to ~ 2Mq and ~ lMg, respectively, and our least 
massive stars (i.e. sink particles) are well above this limit. 
Rather, the effective Jeans mass at T ~ 100K and densities 
of ~ 10 5 — 10 6 cm" 3 prevent the formation of low-mass stars. 

When interpreting our simulation results, there are sev- 
eral caveats that need to be kept in mind. First, our nu- 
merical model does not include shear. Strong shear motions 
may mimic the EOS effects discussed here, as shear adds 
stability and thus requires larger Jeans masses for collapse 
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Figure 3. (a) Mass spectrum of gravitational condensations in the starburst calculation at a time when 15% of the gas is converted into 
collapsed objects (which we identify as direct progenitors of individual stars). To guide the eye, we indicate a slope -1.0 and the Salpeter 
slope -1.3 with dotted lines. The mass function in our simulated starburst environment shows a broad peak in the range 10 Mq — 25 Mq 
and falls off for larger masses. It is thus top-heavy compared to the IMF in the solar neighborhood (Salpeter 1955; Kroupa 2002; Chabrier 
2003). (b) Mass spectrum of collapsed objects in a calculation focusing on nearby molecular clouds (see Jappsen et al. 2005). It agrees 
well with the IMF in the solar vicinity. For comparison we overplot the functional forms proposed by Kroupa (2002) with dashed lines 
and by Chabrier (2003) with dotted lines. Our two calculations differ mainly in the adopted EOS, i.e. in the chemical and thermodynamic 
state of the star forming gas, other parameters are comparable. 



to occur. However, the Arches cluster is bound. Thus the 
Galactic centre shear field cannot play a dominant role in 
the inner parts of the cluster. Second, our numerical model 
does not take the effects of magnetic fields into account, 
which may be of considerable strength in the Galactic cen- 
tre (Yusef-Zadeh & Morris 1987, but also see Roy 2004 for 
lower estimates). However, even if there is a rough equipar- 
tition between kinetic and magnetic energy, the chemical 
and thermodynamic properties of the gas are not strongly 
affected. Our results will still hold at least qualitatively, in 
the sense that an extreme environment leads to deviatiations 
from the standard Salpeter IMF. Third, the use of sink par- 
ticles does not permit us to resolve close binary systems. 
Massive stars in the solar vicinity are almost always mem- 
bers of a binary or higher-order multiple stellar system (e.g. 
Vanbeveren et al. 1998). If this trend holds also for starburst 
environments, then the peak of the stellar IMF will lie below 
the value reported here. For instance, if each unresolved sink 
particle in our calculation separates into a binary star, in a 
statistical sense our mass spectrum needs to be shifted to 
lower masses by a factor of 0.5. Finally, protostellar feedback 
may locally affect the accretion onto individual protostars. 
In this case the mass content of the sink particle may only 
poorly reflect the mass that ends up in a star. However, even 
in the extreme case that half the mass is removed by feed- 
back during collapse (for estimates, see Yorke & Sonnhalter 
2002; Krumholz, McKee, & Klein 2005), deviations from the 
standard IMF will still persist. 

For typical molecular clouds in the Milky less than a 
few percent of their mass takes part in star formation (e.g. 
Myers et al. 1986) and this fraction goes up by a factor of 
a few for cluster-forming cores (e.g. Lada & Lada 2003). 
A number of observations (Paglione et al. 1997; Mooney & 
Solomon 1988) indicate that starburst systems like NGC253 
and M82, and luminous infrared galaxies in general, have a 



larger fraction of their interstellar gas mass at high densities 
(Gao & Solomon 2004). Consequently, their SFE's are up by 
as much as an order of magnitude. Our simulations cover this 
range and the statistics of our mass spectra do not change 
above a SFE ~ 10% in both runs. Hence, the precise SFE 
that pertains to a starburst environment does not influence 
our results as long as it is larger than 10%. 

The computed star formation rate (SFR), defined as the 
change in mass with time of the sink particles, is typically 
860 M Q yr _1 / k P c2 for a SFE > 10 % and when normalised 
to a surface area of 1 kpc 2 , which is roughly the size scale 
of the nuclear region inside a starburst galaxy. This number 
lies well within the fiducial range of 50 - 1000 M Q yr-7kpc 2 
inferred for most starburst systems (e.g. Kennicutt 1998, 
Scoville & Wilson 2004). 

When turning to distant starburst galaxies in the early 
universe, the low-mass cut-off at 7 Mq seen in the simulated 
local starburst region seems at first glance difficult to recon- 
cile with the mass-to-light ratio and the stellar population 
synthesis models inferred from global observations (Kauf- 
mann et al. 2003). However, we emphasise again, that we are 
focusing on an extreme case and on a clearly localised, iso- 
lated region only. In reality these extreme (warm and dusty) 
environmental conditions will not apply to all regions in- 
side a starburst galaxy. There will be pockets of colder gas 
with different (7 < 1) EOS that are less exposed to radi- 
ation (Spaans & Silk 2000) and that behave like Galactic 
star-forming regions. Under these conditions the studies by 
Jappsen et al. (2005) and Larson (2005) indicate that a nor- 
mal, Salpeter-like IMF results. This also suggests that the 
relative contribution of the extreme IMF found in this work 
can be connected directly to the observations. The fraction 
of molecular gas at densities > 10 4 cm -3 that enjoys tem- 
peratures larger than 50 K should be a strong indicator of 
deviations from a Salpeter IMF. Future work will address 
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the issue of stellar population matching and will compare 
our results with observed M/L ratios and warm, high den- 
sity gas mass estimates. 
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